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The statistics of polymers advected by a turbulent flow are investigated. To limit the polymer 
lengths above to coil-stretch transition, a FENE-P type relaxation law is used. The turbulence is 
modeled by a random strain, delta-correlated in time and with Gaussian statistics. The distribution 
of polymer lengths for both the coiled and stretched states are derived, from which are obtained 
analytical expressions for the moments of the distribution. The polymer stress on the fluid decreases 
linearly with the inverse of the Deborah number. The degradation (breaking) of the polymers is 
also discussed, showing that even for large Deborah number some polymers can remain unbroken. 
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One of the principal motivations to study the proper- 
ties of polymers diluted in a turbulent flow is to under- 
stand the Toms effect, also known as drag reduction [1, 2], 
whereby the effect of a small concentration of polymers 
is to dramatically reduce the pressure required to main- 
tain a turbulent flow at a given velocity in a pipe, ft is 
thought that the extension of the polymers by the strains 
in the turbulent flow (the coil-stretch transition) is cru- 
cial to that effect [f, 3]. More recently, there has been 
great progress in understanding the microscopic behavior 
of biological polymers such as DNA [4-6] , which give con- 
fidence that there is a regime where the polymer tension 
dominates the relaxation time [7], as opposed to hydro- 
dynamic effects. 

The dynamics of polymers in incompressible turbu- 
lence was studied recently in Refs. [8—11] . Balkovsky 
et al. [8, 10], building on the approach of Lum- 
ley [12], modeled the turbulence as a homoge- 
neous, isotropic, delta-correlated random strain — the 
Kraichnan-Kazantsev model of turbulence [13-18] (see 
the review by Falkovich et al. [19]) — for polymers 
whose relaxation is described by the Hookean dumbbell 
model [20]. Chertkov [9] extended the analysis to gen- 
eral nonlinear relaxation models. Experiments [11] and 
numerical simulations [21] validate many of the findings 
of these papers. 

In this letter we explore the coiled and stretched states 
of polymers assuming the popular FENE-P model [20, 
22, 23] (Finite Extension Nonlinear Elastic) for the relax- 
ation force of the polymer. The advantage of the FENE 
model is that its nonlinear nature prevents the polymers 
from being infinite in length above the coil-stretch transi- 
tion: it limits them to a finite, fixed length. The results 
for the PDF of the polymer lengths in both the coiled 
and stretched states are presented, as well as moments 
of the distribution, the polymer stress on the fluid, and 
a discussion of polymer degradation by breaking. For 
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definiteness and simplicity, Gaussian statistics for the 
strain are assumed throughout, though the calculation 
herein could be redone for more arbitrary statistics using 
the path integral formalism used by Chertkov [9]. The 
emphasis is thus on presenting closed form, analytically 
simple expressions that give a reasonable description of 
the physics. Compressibility effects are neglected, as is 
the feedback of the polymers on the flow, justified when 
the polymer stress is much less than vX [8], where v is 
the fluid (Newtonian) viscosity and A is the mean rate of 
strain. 



I. FOKKER— PLANCK EQUATION 

The shape of the polymers is characterized by the sym- 
metric conformation tensor, A, which describes the de- 
formation of the polymers from a coiled ball of radius po 
into a stretched ellipsoid in the presence of flow. We 
take the equation of motion for the conformation tensor 
to be that of a general viscoelastic fluid with a relaxation 
force term (only the longest restoring time is taken into 
account), 

A An 

— - A lk d k vP + d k u l A k > - - (f(A kk )A^ - pi 6 ij ) , 

(1) 

where f(A kk ) depends only on the trace A kk of A, and 
is equal to 1 when A = pjjl, and d/dt = d/dt + u k d k is 
the advective derivative. We take the velocity field to be 
incompressible, dkU k = 0. (Unless otherwise noted, re- 
peated indices are summed.) The conformation tensor A 
can be diagonalized, with eigenvalues p\ that evolve ac- 
cording to 

^jp = O-aa Pa~- (f(p0Pf})pa ~ pll Pa) , (2) 

with & a p(t, x) := e l a di~u l e| the velocity gradient tensor 
projected on the orthonormal eigenvectors, e a . There is 
no sum over a on the right-hand side of (2). 

When the viscous scale of the flow is much larger than 
the polymer length (the typical situation), the turbulence 
appears locally as a constant strain. In the context of our 



2 



model, the polymers are stretched ellipsoids whose ma- 
jor axis is aligned with the dominant stretching direction, 
associated with the largest Lyapunov exponent, positive 
here for an incompressible flow. This alignment is ob- 
served in numerical experiments [23]. We thus neglect 
the subdominant directions in (2) and write 

f =Xp-±(f(p 2 )p-pl/p), (3) 

where p is the largest p a and X(t) is the Lyapunov expo- 
nent associated with the dominant direction of stretch- 
ing. 

We take X(t) to be a random variable satisfying 

(X(t)X(t'))-X 2 ^S(t-t')A; (X(t))=X, (4) 

where the angle brackets denote a average over X(t). In 
general we have A < A, since the standard deviation can- 
not be larger than the mean for the positive Gaussian 
variable A. 

To get a Fokker-Planck equation for the PDF of /?, 
we introduce the characteristic function, Z(t; p) — 
(exp(hip)). We can then derive an equation of motion 
for Z from (3) and average. Gaussian integration by 
parts [24-27] allows evaluation of terms of the form (AZ), 
finally yielding a closed equation for Z. Inverse Fourier 
transformation of Z with respect to p then gives the equa- 
tion of motion for 3 3 (t,/j), the PDF of p, 

d t j> = \ a d pP d pP y - x d p P ? + i d p (/( P > - P l/ P ) ?. 

(5) 

(The dp act on all terms to their right.) This is the same 
Fokker-Planck equation as derived by Chertkov [9], ex- 
cept that instead of diffusivity we shall use the minimum 
polymer size po as an ultraviolet cutoff. This also allows 
us to treat the coiled state in Section II. 

II. THE COILED STATE 
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FIG. 1: Moments <p")/po of the PDF of polymer lengths in 
the coiled state for different values of n, with £ := A/ A = 1. 

where j(a, x) is the lower incomplete gamma function. 
(7(a, oo) is equal to r(a), the complete gamma function.) 
For £ > £, the denominator in (7) is positive and finite, 
and for n < 2(£ — £) so is the numerator. 

For n — 2(£ — £) the numerator diverges, so that the 
moment (p n ) becomes infinite. When £ = £ the n = 
moment diverges, which means that the PDF y c becomes 
unnormalizable. Thus the physical picture is as follows: 
for large £ — £ only moments associated with large n di- 
verge, but as £ — £ is lowered more moments diverge, 
indicating that on average more polymers have anoma- 
lously long lengths, that is, they are uncoiled. For £ = £ 
all positive (n > 0) moments are infinite, which indicates 
that the polymers are essentially all in a stretched (un- 
coiled) state, but also pointing to a breakdown in the 
theory because of the infinite moments. 

The combination £ — £ can be rewritten 



We first deal with the case of relatively unstretched 
polymers, so that their effective radius (the maximal ex- 
tension of the ellipsoid describing the deformation of the 
polymers) is much less than the maximal extension of 
the polymer. In that case we assume that the polymers 
are in a linear force response regime and we can use the 
Hookean dumbbell model, f(p 2 ) — 1. Since the determi- 
nant of A is greater than (p^) 2 (see the appendix), we 
can assume that the largest principal axis p (the only 
one we model) is larger that po, and assume the PDF 
of p vanishes for p < po- 

The Fokker-Planck equation (5) with / = 1 can be 
solved for the steady-state distribution, 

5'c(p)^p- 1 - 2( «- C) cxp(-£p2/p 2 ), (6) 

where £ := 1/Ar, £ :— A/A. From the PDF we can 
derive an exact expression for its moments, 

(P n ) 



£-£ = £(De- 1 -l); Do := At, 



(8) 



Po 



cn/2 7 (£-£-n/2,£) 
7(£-C,0 ' 



(7) 



where De is the Deborah number, the ratio of the polymer 
relaxation timescale over the advection timescale. The 
transition to the stretched state thus occurs for De > 1. 
A few moments are plotted CIS Si function of the Deb- 
orah number in Fig. 1, showing their divergence as De 
approaches unity. 

For £ + 1 moderately large (fluctuations are small com- 
pared to the mean), the main contribution to the integral 
representation of the incomplete gamma function occurs 
well below £, so we can replace the upper integration 
bound £ by oo in (7) and obtain 



<£> ^ r /2 r(£ - £ - n/2) 
Pi T(£-£) ' 



(9) 



The even moments then reduce to simple products of 
monomial factors according to the factorial property of 



3 



the gamma function. A simplification occurs for n = 2, 

^^^|^ T = [l-A(A + r)]- 1 . (10) 

This is a particularly significant moment, since the poly- 
mer stress on the fluid is proportional to A ~ p 2 . 

Note that the singular denominator in (10), which 
comes from the gamma function numerator of (9), im- 
plies that the polymer length can be much larger than p 
even below the coil-stretch transition. Of course, when 
it becomes too large we must include nonlinear contribu- 
tions to the relaxation force, the topic of the next section. 



III. THE STRETCHED STATE 

As we saw in Section II, the Hookean dumbbell model 
fails once the Deborah number reaches unity, because 
it cannot prevent the polymers from becoming infinitely 
long [3]. One standard way to prevent this is to use the 
FENE-P model [20, 22, 23], for which 



f(p 2 ) = 



pI 



(11) 



The restoring force of the polymers is now such as to 
limit their length to p m (though they make break before 
reaching that length, as we see below). We assume that 
the elongation of the polymers is smaller than the viscous 
length, so that their length is limited by their own elas- 
ticity and not by the hydrodynamic mechanism of Tabor 
and de Gennes [28, 29]. 

Since in this section we are interested in the stretched 
state of long polymers, we may neglect po *C p m both 
in (5) and (11), which leads to the stationary distribution 



(12) 



Unlike the distribution obtained in [9, 10], this PDF ex- 
hibits a cutoff at p m , reflecting the fact the limited size of 
the polymers. We also impose a cutoff at po for small p, 
because p cannot be less than the equilibrium length p 
(see the appendix). In Rcf. [9] diffusivity was used to 
regularize the PDF at small p, but both methods are 
equivalent since po drops out in our approximation. 
The moments of the PDF (12) are found to be 



(P n ) _ r(C - j + n/2) r(c + On/2 
p« r(c-C)r(c + i + n/2) (C + i) n/2 



(13) 



where (x) n := T(x + n)/T(x) is the Pochhammer sym- 
bol. This expression is valid for ( — £ > 0, or equivalently 
for De > 1. At De = 1, Eq. (13) predicts zero length for 
the polymers, but this is just a symptom of the neglect 
of po. Rather, the coiled solution of Section II must then 
be used. 

For n = 2, Eq. (13) again takes a simplified form, 



n2 



1 - De" 



C + l A + A l + C- 1 



(14) 




FIG. 2: Positive moments {p n )/p™ of the PDF of polymer 
lengths in the stretched state, with £ := A/ A = 1. 




FIG. 3: Negative moments (p") /p™ of the PDF of polymer 
lengths in the stretched state, with £ := A/A = 1. 



so the average polymer stress (proportional to (p 2 )) is 
linear in the inverse Deborah number. Some positive 
(n > 0) and negative (n < 0) moments are plotted as 
a function of the Deborah number in Figs. 2-3. The 
negative moments diverge as De approaches unity, a con- 
sequence of the neglect of p . 

As mentioned in Ref. [9], the PDF of length can 
also be used to estimate the fraction of polymers that 
break. (This is a kind of degradation, though not the 
same as say for DNA, where the two strands become 
separate [30].) Assume a polymer breaks if the ten- 
sion F(p 2 ) = f(p 2 )p/r exceeds a critical value, F c . To 
that critical tension corresponds a critical length, p c , ob- 
tained by solving 



F c = -f(pl) Pc = --^- i p c 

T T Pm - Pc 



(15) 



a quadratic in p c . The fraction of polymers that sur- 
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FIG. 4: Fraction of polymers that survive as a function of 
the inverse Deborah number, for different values of £ :— A/ A 
with pc/pm = 0.7. 
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FIG. 5: Stability diagram showing that D := pip2p3 must 
have value greater than or equal to unity. Constant-D curves 
are dotted, with the D = 1 curve dashed. 



vive is equal to the fraction shorter than p c , obtained by 
integrating the PDF (12) from to p c , 



Prob(p < p c ) = 



s(c-£,e+i) ' 



(16) 



where B(a,b) :=T(a)T(b)/T(a + b) is the beta func- 
tion, and B(z; a, b) is the incomplete beta function. 
Since 5(1; a, 6) = B(a,b), for p c = p m none of the poly- 
mers break. The fraction of surviving polymers (16) is 
plotted in Fig. 4 for different values of (. For large £ 
(small fluctations), the survival probability becomes a 
step function with "knee" at De _1 = 1 — (p c /pm) 2 : ei- 
ther all polymers survive or all are broken. 

A consequence of (16) is that a significant frac- 
tion of polymers can survive even for large Deb- 
orah number; For ( C ( (Dc 3> 1) and ^ 1 , 
Prob(p < p c ) ~ (pc/pm) 2 '', which is not small if £ is not 
too large. In other words, some polymers survive break- 
age for large A if the fluctuations in A (given by A) are 
also large. Of course, this treatment is for one correla- 
tion time of the turbulence. In reality a given polymer 
is likely to be exposed to many different random strains. 
Because we have assumed the turbulence to be delta- 
correlated, the probability of survival of the polymers de- 
cays roughly as (p c /pm) 2< ' N , where N = t/t COIT with t the 
elapsed time and t corr the correlation time. (More gen- 
erally, for smaller Deborah number, take the Nth power 
of (16).) 



so we let D := pxpiPz- (We rescale p such that ||po|j = 1-) 
For small volumes (low stretching) the Hookean approxi- 
mation applies, so we let / = 1 in (2); the quantity logD 
obeys 



r-log2) = -3+( i0 r ; 



P2 2 



-pf)=:L(p). (17) 



The sign changes of the derivative of log D are the same 
as those of D since the logarithm is monotonically in- 
creasing. Figure 5 shows the sign of £ in the pi—ps 
plane, assuming pi = 1. Without loss of generality we 
also assume pi > ps, so we need not consider the shaded 
area. The solid curve is £ = 0, where the determinant 
is constant. The dotted curves show lines of constant D, 
the dashed line being D = 1. The dotted and dashed 
curves meet at the spherical solution. Below the solid 
curve, T> must increase, and above it must decrease. The 
line D = 1 lies entirely below L — and is tangent to it 
at pi = P3 = 1. Assuming that the polymers are initially 
in the spherical (rest) state, we see that D must either re- 
main at or increase away from unity. In fact the only dis- 
turbances that do not linearly decay back to the spherical 
state are those along the solid line. The picture for pi ^ 1 
is three-dimensional but is qualitatively the same. We 
conclude that D > 1, so that the volume of the ellipsoid 
is no smaller than the equilibrium sphere. Hence, the 
largest polymer dimension, pi, is never smaller than po, 
justifying its use as a small-pi cutoff. 



APPENDIX: LOWER BOUND ON POLYMER 
VOLUME 
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